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ABSTRACT 

We present and analyse optical and ultra-violet (UV) colours for passive and optically-red 
Coma cluster galaxies for which we have spectroscopic age and element abundance estimates. 
Our sample of 150 objects covers a wide range in mass, from giant ellipticals down to the 
bright end of the dwarf-galaxy regime. Galaxies with ongoing star-formation have been re- 
moved using strict Ha emission line criteria. We focus on the colours FUV — i, NUV — i, 
FUV — NUV, u* — g and g — i. We find that all of these colours are correlated with both 
luminosity and velocity dispersion at the > 5<r level, with FUV — i and FUV — NUV 
becoming bluer with increasing 'mass' while the other colours become redder. We perform 
a purely empirical analysis to assess what fraction of the variation in each colour can be ac- 
counted for by variations in the average stellar populations, as traced by the optical spectra. 
For the optical colours, u* — g and g — i, most of the observed scatter (^80 per cent after 
allowing for measurement errors and for systematic errors in u* — g) is attributable to stel- 
lar population variations, with colours becoming redder with increasing age and metallicity 
(Mg/H). The FUV — i colour becomes bluer with increasing age and with increasing Mg/H, 
favouring the 'metal-rich single-star' origin for the UV upturn. However, correlations with 
the optically-dominant stellar populations account for only about half of the large observed 
scatter. We propose that the excess scatter in FUV — i may be due to a varying proportion of 
ancient stars in galaxies with younger (SSP-equivalent) average ages. The NUV — i colour 
is sensitive to SSP-equivalent age and Mg/H (in the same sense as optical colours), but also 
exhibits excess scatter that can be attributed to 'leakage' of the FUV-dominant old hot popu- 
lation. After applying a correction based on the FUV — i colour, the much of the remaining 
variance in NUV — i is attributable to variations in the spectroscopic parameters, similar to 
the results for optical colours. Finally, the FUV — NUV colour is surprisingly well behaved, 
showing strong correlations with age and metallicity, and little residual scatter. Interpreting 
this colour is complicated however, since it mixes the effects of the main-sequence turn-off, 
in the NUV, with the variation in the hot post-RGB content dominating the FUV. 

Key words: galaxies: elliptical and lenticular, cD - galaxies: stellar content - ultraviolet: 
galaxies 



1 INTRODUCTION 

The optical and near-infrared luminosities of old ( > 3 Gyr) stellar 
populations are dominated by main sequence, sub-giant and red 
giant branch (RGB) stars, with a combined contribution of ~80 per 
cent (Maraston 2005). In the ultraviolet (UV) regime, by contrast, 
these fairly well-constrained evolutionary phases are faint, and the 
total light output becomes sensitive to hot post-RGB phases that 
are much less well understood. 

In particular, the spectral energy distributions of many giant 
ellipticals and spiral bulges exhibit an 'upturn' in the far-UV (FUV) 
below ~2000A (Code & Welch 1979). Early UV spectroscopy 
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(e.g. Welch 1982, Ferguson et al. 1991) and imaging (e.g. Bohlin 
et al. 1985; O'Connell et al. 1992) showed that the upturn was 
not due to young massive stars, and hence low-mass evolved stars 
in old stellar populations are required (see O'Connell et al. 1999 
for an extensive review). The favoured candidate sources for the 
UV upturn are extreme horizontal branch (EHB) stars, and their 
hot post-horizontal-branch descendants (Greggio & Renzini 1990), 
which have temperatures T c g~25 000K. Such stars are observed 
in globular clusters, and also in the field, where they are desig- 
nated sub-dwarf B stars (see Catelan et al. 2009 for a review of the 
EHB among other horizontal branch phenomena). Their formation 
is thought to require either large mass loss on the RGB, occurring 
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at high metallicitj0 (Greggio & Renzini 1990; Yi, Demarque & 
Oemler 1997) or through a variety of binary interactions (Han et al. 
2007). An unresolved question is the degree to which these rather 
exotic stellar phenomena are related, on a broad scale, to the prop- 
erties of the stellar populations which dominate the optical light of 
the same galaxies. 

Observationally, one way to explore the systematic variation 
in FUV output for (optically) red galaxies is to compare UV colours 
against luminosity and velocity dispersion, in analogy with the tight 
scaling relations obtained for optical colours. However, a more di- 
rect test for the physical origin of the UV flux is to compare against 
indicators for the age and metallicity of the stellar populations that 
dominate in the optical. Optical colours and individual absorption 
line-strengths (e.g. Mg6 or H/3) provide some constraints on the 
stellar populations, but are always sensitive to a combination of 
both age and metallicity effects. An alternative method, pursued in 
this paper, is to correlate the UV colours against ages, metallicities 
and element abundance ratios derived from line indices, through 
comparison to spectral synthesis models. This is by no means the 
first attempt to address these questions, and we review some rele- 
vant previous work in the following paragraphs. 

Using observations of 31 early-type galaxies with the Inter- 
national Ultraviolet Explorer (WE), Burstein et al. (1988) showed 
that FUV — V colour becomes bluer with increasing velocity dis- 
persion, a, as the upturn component strengthens in the most mas- 
sive galaxies. Moreover, the same colour was found to be anti- 
correlated with the Lick Mg2 absorption index. Since this index 
is primarily a metallicity indicator (although also sensitive to age) 
Burstein et al. concluded that metallicity is the fundamental pa- 
rameter underlying the observed correlations. O'Connell (1999) 
repeated this test with a larger set of line indices from Trager et 
al. (1998), and found that FUV — V shows clear trends with in- 
dices tracing light-element abundances, but is essentially uncorre- 
cted with indices dominated by iron abundance. 

A dramatic improvement in the systematic study of UV emis- 
sion from passive galaxies was enabled by the Galaxy Evolution 
Explorer (GALEX) satellite (Martin et al. 2005; Morrissey et al. 
2007), with its wide field of view and two-channel imaging capabil- 
ity (NUV: A cff = 2310 A FUV: A cff = 1530 A). In an early study 
of UV colours of 'quiescent early-type' galaxies with GALEX, Rich 
et al. (2005) concluded that there was 'no correlation between [the 
FUV — r colour] and any parameter sensitive to the global metal- 
licity', including the Mg2 index used by Burstein et al. (1988). The 
Rich et al. sample was limited to luminous galaxies, with essen- 
tially no coverage at a < 100 km s _1 . Other early GALEX studies 
by Boselli et al. (2005) and Donas et al. (2007) focused on nearby 
galaxies, and covered a wider range in luminosity. Both works re- 
covered correlations between FUV — NUV and Mg2, in the sense 
observed by Burstein et al. (1988), although Donas et al. show that 
the correlation is much weaker for FUV — V. They emphasise that 
the tight correlation in FUV — NUV is likely driven by metal- 
line blanketing in the NUV from the main-sequence turn-off stars, 
rather by than variation in the hot UV-uptum sources. 

More recently, Bureau et al. (2011, hereafter Bll) have re- 
visited this issue using GALEX colours combined with SAURON 

The formation of EHB stars in low-metallicity populations appeared un- 
able to reproduce the upturn without adopting unrealistically high ages. 
Chung et al. (2011) have recently argued that helium-enhanced sub- 
populations (for which there is independent evidence in some globular clus- 
ters) can develop a strong upturn at ages less than the Hubble time, since 
main sequence lifetimes are shorter for higher helium abundances. 



integral-field spectroscopy, for a sample of 48 very nearby early- 
type galaxies. Their analysis has an advantage over previous work, 
in that colours and spectroscopic parameters can be extracted at 
identical apertures. After excluding galaxies with strong H/3 ab- 
sorption, indicative of recent star-formation, Bll find that galaxies 
with bluer FUV — V have stronger Mg6, similar to the original 
result of Burstein et al. (Note that Bll follow Burstein et al. in 
treating colour as the independent variable in the regression). They 
also find that galaxies with bluer FUV — V have lower H/3 absorp- 
tion, suggesting older populations or higher metallicity, while there 
is no clear trend for Fe5015, an iron-dominated index. They con- 
clude that Q-element abundance or total metallicity is an important 
driver of the UV upturn. 

Although Mgfe (like Mg2) is considered primarily a metal- 
licity indicator, it also depends on age. Hence a correlation of 
FUV— optical colour with this index does not unambiguously sig- 
nal that the UV upturn depends on metallicity. Comparison of the 
UV colours against stellar population parameters (age, metallicity, 
etc) estimated from multiple indices should help to separate the 
effects of the physically meaningful quantities. This was attempted 
for the NUV — J colour by Rawle et al. (2008), who concluded that 
metallicity effects were dominant in driving the systematic vari- 
ation in this colour, with a residual dependence on age, while a 
large (0.25 mag) excess scatter remained unexplained by the cor- 
relations. Loubser & Sanchez-Blazquez (2011) performed similar 
comparisons for a sample of brightest cluster galaxies (BCGs) and 
a control sample of similarly massive ellipticals that are not BCGs. 
They found no significant correlations of FUV — NUV either with 
linestrength indices (including Mgb5177), or with age, metallicity 
(Z/H) or abundance ratio a/Fe taken individually. The construction 
of their sample, by its nature, limits the baseline in mass and stel- 
lar population properties which would be necessary to constrain 
any such correlations. Carter et al. (2011) have applied a similar 
method for a sample of nearby ellipticals with spectroscopic pa- 
rameters compiled from a range of literature sources. Testing for 
correlation with each single parameter in turn, they found signifi- 
cant negative correlations of FUV — NUV with total metallicity 
(Z/H), and with a/Fe but not with iron abundance (Fe/H); the age 
dependence was not explored. 

In this paper, we use deep GALEX imaging of the Coma clus- 
ter, together with optical data from the Canada-France-Hawaii 
Telescope (CFHT) and the Sloan Digital Sky Survey (SDSS), 
to measure ultraviolet and optical colours for a sample of red- 
sequence galaxies for which ages, metallicities and abundance ra- 
tios have been determined from high-S/N spectroscopy. The key 
improvements over previous works are (1) the use of (in principle) 
meaningful stellar population parameters, rather than individual in- 
dices which mix age and metallicity effects, and (2) application 
to a sample spanning a wide range in galaxy mass (a = 50 — 
400 km s _1 ). We analyse the correlations of the colours with spec- 
troscopic ages and abundances, making use of multi-parametric fits 
to disentangle the effects of each parameter. In particular, we aim 
to determine what fraction of the variation in each colour is driven 
by variations in the stellar populations, as constrained from the op- 
tical, and what fraction is 'excess scatter', apparently unrelated to 
the optically-dominant populations. Our approach is purely empir- 
ical, and independent of any population synthesis models for the 
UV colours. 

The structure of the paper is as follows: the optical spec- 
troscopy and the GALEX and optical imaging data are described 
in Section [2] The results from fitting the colours with single- and 
multi-parameter models are presented in Section [3] together with 
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an exploration of residual trends and tests for effects of systematic 
errors. We discuss the implications of our results in Section [4] fo- 
cusing especially on the questions of the origin of the UV upturn, 
and whether the scatter in the NUV implies widespread recent star 
formation among optically passive galaxies. Our conclusions are 
summarised in Section[5] 



2 DATA 

Our approach is to measure UV/optical colours for galaxies drawn 
from a catalogue of galaxies with measured spectroscopic ages and 
metallicities. Full details of the spectroscopy are given in related 
papers (Smith et al. 2012; Smith et al., in preparation). Broadly 
the parent sample comprises spectra for 242 bright Coma cluster 
members from SDSS (analysed in Price et al. 2010) and 169 fainter 
members observed using long integrations with Hectospec at the 
6.5m MMT. The MMT sample is an extended version of that re- 
ported by Smith et al. (2009). Velocity dispersions were compiled 
from a variety of literature sources, supplemented with new ob- 
servations from VLT/FLAMES, and combined using observations 
from multiple data-sources to determine relative systematic offsets. 

The combined MMT and SDSS sample has been homoge- 
neously re-analysed, to measure absorption line strength indices 
and emission line equivalent widths (used in defining a passive 
galaxy sample). Ages, metallicities (Fe/H), and abundance ratios 
(Mg/Fe, Ca/Fe, C/Fe and N/Fe) were measured via comparison to 
the Schiavon (2007) simple stellar population (SSP) models, using 
a new model inversion code. The first step in this process is a non- 
linear optimisation to derive Fe/H and age assuming solar abun- 
dance ratios, using the observed H/3 and (Fe) = i(Fe5270+Fe5335) 
indices. Then we use a fast linear fit (i.e. a matrix inversion) to 
estimate a correction to the abundance pattern required to match 
additional indices (Mg6, Ca4227, C 2 4668 and CN 2 ), at the fitted 
age and Fe/H. Using this updated abundance pattern, we repeat the 
non-linear fit to obtain an improved estimate of the age and Fe/H. 
Alternating between these steps, the process converges within a 
few iterations. Parameter errors (and covariances) are obtained by 
propagating the index errors through the fitting process using mon- 
tecarlo simulations (50 realisations per galaxy). We allow extrap- 
olation of predictions beyond the model limits (e.g. to unphysi- 
cally large ages), in order to characterise the errors correctly. We 
have confirmed that our method yields results consistent with the 
much slower method of Graves & Schiavon (2008), when applied 
to the same data. The derived properties (especially age) should 
be interpreted as 'SSP-equivalent parameters', i.e. they are the pa- 
rameters of the single-burst, single-metallicity, single-abundance- 
mixture population which best reproduces the measured indices. 

We draw the UV data from proprietary and archival GALEX 
observations in NUV and FUV There are two deep observations: 
GI5_025001_COMA with 15ksec integration in the cluster core 
(see also Smith et al. 2010) and GI2_046001.COMA3 with 30ksec 
exposure in a field centred 0.9 degrees to the south-west (see also 
Hammer et al. 2010). Retrieving additional archival tiles with inte- 
gration times 1-5 ksec from the archive resulted in coverage for al- 
most all galaxies in our spectroscopic sample. For galaxies covered 
by more than one GALEX tile, we employ only the deepest avail- 
able observation. Deep observations with GALEX are confusion- 
limited, and many of the faintest objects in our sample would not be 
individually detected by running, e.g. SExtractor (Bertin & Arnouts 
1996) on the GALEX images. Instead, we simply extract the UV 
fluxes from circular apertures with centres fixed at the optical coor- 



dinates. We adopt an aperture of lOarcsec diameter, corresponding 
to ~5kpc at the distance of Coma, which is sufficiently large to 
avoid sensitivity to variation in the GALEX point-spread function 
(5-6 arcsec FWHM). A local sky estimate is obtained from an an- 
nulus of radius 30-60 arcsec, with contaminating sources identified 
by k-sigma clipping the (deeper) NUV images, and subsequently 
excised from the sky statistics in both bands. 

For optical photometry in redder passbands, we use aperture 
magnitudes from the 7th Data Release of the Sloan Digital Sky Sur- 
vey, for the g and i bands, interpolating from the published aper- 
tures to the 10 arcsec diameter used for the UV. Since the SDSS Ri- 
band photometry is very shallow, we instead use dedicated observa- 
tions obtained with MegaCam at the Canada-France-Hawaii Tele- 
scope (CFHT). The CFHT u*-band imaging covers a 2.8x2.8 deg 2 
area with a grid of nine slightly-overlapping pointings. The u* filter 
is ~200 A redder than SDSS u. The individual exposures for each 
~1 deg 2 pointing were processed and stacked at the observatory us- 
ing the Elixir pipeline (Magnier & Cuillandre 2004). Magnitudes in 
10-arcsec apertures were measured from the stacks, which have a 
total exposure time of 1360 sec. 

The adopted random errors in the aperture magnitudes and 
colours include contributions from photon statistics and from con- 
fusion noise. We estimate the latter in all five bands, by measuring 
the 68 percentile interval of the total flux in 10-arcsec apertures 
placed randomly in a representative image, and propagating the ef- 
fects into errors as a function of galaxy magnitude. Confusion noise 
is the dominant source of random error in the UV for the deepest 
GALEX tiles, due to the low spatial resolution of the images, and 
also because the cluster galaxies are relatively faint in the UV com- 
pared to the background sources. Even in the optical, however, con- 
fusion can be comparable to the (generally very small) photon er- 
rors, given the large apertures used for the photometry. The adopted 
errors do not include contributions from systematic calibration er- 
ror, since many of the galaxies in the sample are drawn from only 
a small number of observations (i.e. GALEX tiles, SDSS drift-scan 
runs, MegaCam pointings) per data source, and hence calibration 
errors affect large parts of the sample coherently. We assess the im- 
pact of systematic calibration errors in Section [3.4.1l bv explicitly 
allowing additional 'offset' terms in our analysis. 

No corrections for galactic extinction are applied. At the high 
galactic latitude of Coma, the mean reddening from Schlegel et al. 
(1998) is (E(B — V))=0.01 mag, with a galaxy-to-galaxy scatter 
of 0.001 mag. Even in the NUV (the band most sensitive to extinc- 
tion), this would introduce a scatter of < 0.01 mag in flux, which is 
negligible in comparison to the observed scatter (~0.2 mag). No k- 
corrections are applied, since they are necessarily model-dependent 
and very uncertain in the UV bands. Although all sample galaxies 
lie at approximately the same distance, their recession velocities 
span a range 3800-9400 kms -1 . Hence band shifting effects may 
be relevant, for example in the u* band, where the long- wavelength 
filter cut-off is close to the steep 4000 A break. We address this 
point in Section 13.4.21 by testing fits which explicitly include a 
redshift-dependent term. 

From the full spectroscopic dataset, we restrict the sample 
to galaxies having age greater than 2Gyr, error less than 0.3 in 
logtssp and [Mg/H], and FUV error less than 0.3 mag. Further- 
more, we retain only galaxies with measured velocity dispersions, 
and with <0.5 A equivalent width in Ha emission (after separation 
from the absorption line). The purpose of the Ha cut is primarily 
to remove galaxies in which the stellar H/3 absorption line, used in 
the estimation of ages and metallicities, is contaminated by nebular 
emission. In addition, the selection against Ha emission helps to 



© 2011 RAS, MNRAS 000.fTMT6l 



4 Russell J. Smith et al. 



remove objects whose UV colours are influenced by ongoing star 
formation in their central regions. (We address the possibility of 
star formation beyond the spectroscopic fibre in Section [3.4.3l ) We 
automatically reject galaxies with UV flux centroids offset from the 
nominal galaxy position by more than 1.5 arcsec, to select against 
objects contaminated by neighbouring sources. (The centroid shift 
threshold is supported by simulations of isolated poisson-sampled 
PSFs, tuned to the depth of the shallower tiles). Finally, all remain- 
ing sources were inspected visually on the GALEX and optical im- 
ages, and five galaxies manually removed from the sample due to 
blending undetected by the centroid offset, or location within the 
halo of a UV-bright star. 

The resulting sample comprises 150 galaxies. For 90 objects 
(60 per cent), the UV measurements are from the deep Coma Core 
tile, while for 25 (17 per cent), the UV data are from the deep 
Coma South West tile. The remaining 35 objects (23 per cent) are 
from the shallower tiles. The merged sample is mainly (97 per cent) 
comprised of galaxies with SDSS spectroscopy, since the UV flux 
errors are typically too large for the fainter galaxies observed that 
were observed with MMT The final dataset used for the analysis in 
this paper is reported in Tables QJ[2] Note that the derived parame- 
ters for SDSS galaxies differ slightly from those tabulated in Price 
et al. (2010) due to different conventions and methods adopted in 
reanalysing the index data. 



determined by minimizing the scatter in colour, this approach dif- 
fers from some previous analyses (e.g. Bl 1). 

We quantify the fits primarily through the 'coefficient of deter- 
mination', R 2 , defined as the fraction of total variance in the colour 
which is accounted for by correlations with the predictor variable: 

_ ^ _ ^2j(yi ~ Pi) _ ^ _ ^residual 
T,i(Vi-V) 2 ^total 

Here, g/j are the measured colours and pi are the colours predicted 
from a linear model, e.g. in the case of a fit to SSP-equivalent age. 
As indicated in the second expression, 1 — R 2 is the ratio of the 
residual variance around the fit to the total variance in the colour. 

We also use a generalisation of the above to account for mea- 
surement errors: 

E 2 2 / 2 2 \ 

s _ { j i e y,i e p,i 1 — | °' rcsidual _ '-'errors \ 

T,i(Vi-V) 2 V ^otal °tot a l / 

where e Vt i are the errors in measured colour, and e Pt i are the errors 
in the predicted colour, e.g. e p = atei ogtsS p in the case of the 
fit to age. T 2 thus represents the fraction of the total variance (in 
observed colour) which is attributable to a combination of (i) the 
correlations with the predictor variables and (ii) the measurement 
errors, both in colour and in the predictors. 



3 LINEAR ANALYSIS OF THE COLOURS 
3.1 Methodology 

In this paper we take a purely empirical and statistical approach to 
assessing the extent to which variations in various physical param- 
eters can account for the variation in the UV and optical colours of 
non-star-forming galaxies. We make no use of population synthesis 
models to predict the colours. Instead, we perform linear regression 
analyses, seeking to find the predictor variable (or variables) which 
best reproduce the observed colours. 

We restrict attention to a subset of colours chosen to highlight 
the behaviour in the UV in comparison to the optical, with little 
redundancy: 

(i) FUV — i : directly probes the UV-upturn sources, measuring 
the flux from old hot stars versus that from old cool stars. 

(ii) NUV — i : the NUV flux likely has contributions from the 
MS turn-off as well as from the upturn sources. 

(iii) FUV — NUV : a colour measurable from GALEX alone. 

(iv) w* — g : a blue-optical colour, measuring the 4000 A break, 
to test the transition between UV and optical regimes. 

(v) g — i : a broad optical colour for comparison. 

We fit the observed colours to the following linear models in 
which the predicted colour is p, given by: 

p = a , L + cll log(L r /10 10 L rjQ ) 

p = fflo.o- + a a log(er/100 kms -1 ) 

P = ao.t + at log(T S sp/10 Gyr) 

P = o ,z + az[Mg/H] 

P = ao,tz + at log(Tssp/10 Gyr) + az[Mg/H] , 

where a is velocity dispersion, L r is total r-band luminosity, Tssp 
is spectroscopic age and [Mg/H] is spectroscopic Mg abundance, 
which is close to the total metallicity, Z/H. The coefficients a are 



3.2 Observed colours versus luminosity and velocity 
dispersion 

We begin by briefly considering the colour-luminosity and colour- 
a relations, as shown in FigureQ] The slopes and other parameters 
of the fits are reported in Table[3] 

The colour-luminosity relations (upper panels of Figure Q} 
show the familiar results that FUV - i and FUV - NUV be- 
come bluer with increasing luminosity, contrary to the behaviour of 
colours at longer wavelengths. Galaxies with FUV — NUV < 0.9 
have a 'classical' UV upturn, in the sense of having rising spectra 
(/a-vs-A) in the UV regime. Such colours are seen mainly at the 
bright end of the luminosity distribution. Our approach is to treat 
the galaxy population continuously, rather than assign any particu- 
lar significance to this value. For reference however, we distinguish 
the 14 objects with FUV — NUV < 0.9 using different symbols 
in this and subsequent figurefl 

The colour-luminosity relations show that the UV colours ex- 
hibit much greater scatter (0.2-0.3 mag) than the optical colours 
(~0.05 mag). In part, this difference in the absolute value of the 
scatter is trivially accounted for by the much larger total colour 
range spanned by the sample. The quantities R 2 and T 2 scale this 
out by comparing the variance around the fit to the total variance 
in colour. The R 2 values indicate that the luminosity dependence 
accounts for only ~15 per cent of the total variance in the 'UV-to- 
optical' colours FUV — i and NUV — i. Luminosity correlations 
account for a larger fraction of variance in the 'optical-to-optical' 
colours u* —g (25 per cent) and g~ i (36 per cent). Perhaps surpris- 
ingly, the luminosity trend also accounts for 30 per cent of the scat- 
ter in the 'UV-to-UV colour, FUV - NUV. We will find through- 
out this paper that FUV — NUV follows tighter correlations than 
the UV- to-optical colours. The T 2 values show that measurement 



2 We note that Yi et al. (201 1) found no bright Coma ellipticals satisfying 
FUV — NUV < 0.9. The apparent discrepancy is probably due to the 
^-corrections applied in that work but not here. 
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Table 1. The galaxy properties used for the analysis in Section [5] The identifications are from Godwin, Metcalfe & Peach (1983), where available, with 
coordinates derived from matching to SDSS. The parameters used as predictors for the colours are the r-band luminosity from SDSS (L r ), the velocity 
dispersion (a) from the compilation in Smith et al. (2012), and the SSP-equivalent age and metallicity (Tggp and [Mg/H]) from Smith et al. (2012). The 
errors on the latter parameters are highly anti-correlated. The error coefficient p crr is needed when computing errors on combinations of age and metallicity. 
The complete contents of this table are provided in the online version of the journal. 



Galaxy R.A. Dec log(L r /L©) log(o-/km s" 1 ) log(T SS p/10 Gyr) [Mg/H] p, 
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Table 2. The 10-arcsec diameter aperture magnitudes from which colours were derived for the analysis in Section[3] The FUV and NUV magnitudes are 
from GALEX, g and i are from SDSS (interpolated from DR7 aperture photometry), and u* is from CFHT (this passband is centred ~200 A redder than SDSS 
u). All magnitudes are on the AB system. The quoted errors include a contribution from confusion noise, but not for calibration uncertainty. Neither extinction 
corrections nor fc-corrections have been applied. The complete contents of this table are provided in the online version of the journal. 
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Figure 1. UV and optical colours as a function of optical luminosity L r (upper panels) and velocity dispersion cr (lower panels). The solid line in each panel 
is the linear fit as recorder in Table[5] The header line indicates the rms scatter about the fit, the coefficient of determination R? (the fraction of total variance 
in the colour that is accounted for by the correlation), and the T 2 parameter which is equivalent to R? but includes the effect of measurement errors. Galaxies 
marked with triangular symbols have FUV — NUV < 0.9, indicating that they are classical UV-upturn galaxies. 



errors contribute an additional ~ 10 per cent of the total variance in 
the UV colours, but make negligible impact on the optical scatter. 

As noted in previous work (e.g. Marino et al. 2011; Carter 
et al. 201 1), the trends of UV colours with velocity dispersion are 
tighter than those with luminosity. The scalings with velocity dis- 
persion, combined with the measurement errors, account for around 
70 per cent of the variance in g — i and FUV — NUV, but only 
30-50 per cent of the variance in FUV — i, NUV — i and u* — g. 
Fitting simultaneously to a and L r (not given in the table or figures) 
indicates there is no significant dependence on luminosity at fixed 
velocity dispersion, consistent with results in the optical (Bernardi 
et al. 2005; Graves et al. 2009). 

To summarize, all colours are strongly correlated (at the > 5a 
level) with optical luminosity and velocity dispersion. These cor- 
relations account for a much of the observed variance in the g — i 
and FUV — NUV colours, but significant excess scatter remains 
in u* — g and in the UV-to-optical colours. 

3.3 Observed colours versus age and metallicity 

We assume that the trends of colour with luminosity and velocity 
dispersion result from underlying variation in the average stellar 
populations in galaxies as a function of these quantities. In this sec- 
tion we use the spectroscopic ages and metallicities to explore the 
degree to which such correlations can explain the observed colour 
variation. 

In the absence of photometric errors (and of internal reddening 
or non-stellar flux sources such as active nuclei), the colours should 
be completely determined from the stellar populations present. In 
practice however, we do not have full information about these pop- 
ulations. The spectroscopic data provide only a single characteristic 
age and metallicity, and estimates for abundance ratios of a few el- 
ements. At best, these are only moments of the true distribution in 
these parameters (e.g. Tssp is weighted to the age of the most re- 
cent star-formation episodes). An important caveat, in the context 
of this paper, is that if the UV colours are driven by sub-populations 
far from the average population, such as the very oldest or youngest 



stars, or the low- or high-metallicity tail' of the metallicity distri- 
bution, then this may not be reflected cleanly in correlations with 
average SSP-equivalent properties. However, spectroscopic age and 
metallicity are still more likely to reflect meaningful information 
about the galaxies than using individual line index strengths, which 
always depend on combinations of these parameters, as well as 
on individual element abundances. At worst, the SSP-equivalent 
parameters may be influenced by effects that are not adequately 
taken into account in the analysis used to derive them. For exam- 
ple, the measured 'age' may partly depend on the incidence of blue 
horizontal-branch stars (e.g. Lee, Yoon & Lee 2000; Schiavon et 
al. 2004; Percival & Salaris 2011), as well as the main-sequence 
turnoff mass. Even in this case, however, the optical age and metal- 
licity provide a useful compression of the spectroscopic constraints 
on the stellar content that are available from the optical light. 

To test for correlations between the colours and the stellar pop- 
ulations which dominate the spectra, we fit colour versus age and 
metallicity (Mg/H), first individually, and then in linear combina- 
tion. We use Mg/H as the primary metallicity tracer, rather than 
Fe/H, since the abundance of a elements better reflects the total 
metallicity Z/H (due to the large contribution of oxygen to the total 
number fraction in metals). The results are shown in Figure [2] and 
the fit results summarised in Table[3] 

The first row of panels in Figure [2] show that it* — g, g — i 
and NUV — i all become redder with increasing age, while 
FUV — NUV and FUV — i become bluer with increasing age. 
The coefficient of age is significant at the >6<r level for all colours 
except NUV — i, where the correlation is only marginally signif- 
icant (2.4<t). In FUV — i, we note an apparent steepening of the 
age dependence among the oldest galaxies. For T ssp >9Gyr (the 
median age of the sample galaxies), the slope is —1.50 ± 0.35 mag 
per decade in age, while for T ssp <9 Gyr the slope is shallower and 
barely significant, at —0.39 ± 0.21 mag per decade. According to 
R 2 (i.e. without accounting for measurement errors), the linear age 
trends make 20-40 per cent contributions to the total variance in 
the colours (except NUV — i). However, the error contribution is 
substantial, because age has fairly large uncertainty, compared to 
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Figure 2. UV and optical colours as a function of age (Tggp) and metallicity (Mg/H). Annotations are as in Figure[T] In the third row, the colours are plotted 
against the optimal linear combination of age and metallicity, as determined from the simultaneous fits. 
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Figure 3. Maps of the average galaxy colour as a function of age and metallicity, created with a Kriging model (e.g. Cressie 1993). Black points indicate 
the galaxies used to generate the maps. The contours are located at deciles of the colour distribution. The upper-right shows the results for NUV — i after 
correcting for UV-upturn 'leakage' in the NUV band (see Section l-OV 
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Table 3. Parameters of linear fits to the observed colours. For each colour, the six lines give the results for (1) a null model fitting only a mean value; 
(2) a fit against luminosity; (3) a fit against velocity dispersion; (4) a fit against spectroscopic age; (5) a fit against spectroscopic metallicity; and (6) a 
simultaneous fit to age and metallicity. The coefficients a refer to a general model of the formp = ao + a^ log(L r / (10 10 Lq)) + a CT log(cr/100 km s -1 ) + 
at log(tgsp / 10 Gyr) + az [Mg/H] , although only subsets of these terms are used in each fit. For each line, the quantity R 2 is the ratio of the total variance 
explained by the fitted model. T 2 is the fraction of variance explained by the fitted model and the nominal measurement errors in both the colours and the 
predictor variables. 
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— 


0.153 


0.20 


0.28 


6.420 ±0.015 






±0.158 ±0.067 




0.168 


0.04 


0.14 


6.272 ±0.031 








+0.612 + 0.126 


0.159 


0.14 


0.38 


6.271 ±0.030 


— 


— 


+0.221 ±0.062 


+0.691 + 0.124 


0.153 


0.20 


0.37 


NUV 1.275 ±0.012 





— 


— 


— 


0.382 


— 


— 


1.287 ± 0.026 


-0.664 ±0.083 


— 


— 


— 


0.320 


0.30 


0.42 


1.613 ±0.026 


— 


-2.177 ± 0.123 


— 


— 


0.216 


0.68 


0.81 


1.185 ±0.026 


— 


— 


-1.145 ± 0.121 


— 


0.302 


0.38 


0.76 


1.578 ±0.069 








-1.364 ± 0.282 


0.355 


0.14 


0.42 


1.580 ± 0.046 






-1.312 ±0.097 


-1.835 + 0.192 


0.237 


0.62 


0.83 


1.436 ±0.019 










0.077 






1.434 ± 0.005 


+0.124 ±0.017 








0.067 


0.26 


0.26 


1.382 ±0.007 




±0.353 ± 0.033 






0.058 


0.44 


0.45 


1.449 ±0.006 






±0.168 ±0.028 




0.069 


0.20 


0.34 


1.370 ±0.014 








±0.300 ± 0.056 


0.071 


0.16 


0.36 


1.369 ±0.011 






±0.202 ± 0.024 


±0.373 ± 0.047 


0.058 


0.44 


0.53 


1.183 ±0.005 










0.048 






1.181 ±0.003 


+0.091 ±0.010 








0.038 


0.36 


0.38 


1.141 ±0.003 




±0.273 ±0.015 






0.026 


0.69 


0.71 


1.194 ±0.003 






+0.137 + 0.015 




0.038 


0.35 


0.61 


1.140 ±0.008 








+0.191 ±0.034 


0.043 


0.17 


0.40 


1.140 ±0.006 






±0.160 ±0.012 


±0.249 ± 0.023 


0.029 


0.63 


0.74 



Table 4. Comparision of fits using age and Mg/H to those using age and Fe/H. An intercept term has also been fitted but is not reported here. 



Fitting age and Mg/H Fitting age and Fe/H 



Colour 


at 


az 


rms 


R 2 


at 


az 


rms 


R 2 


FUV - i 


-1.092 ±0.119 


-1.144 ±0.236 


0.291 


0.38 


-1.092 ±0.159 


-0.267 ±0.250 


0.312 


0.29 


NUV - i 


+0.221 + 0.062 


±0.691 ± 0.124 


0.153 


0.19 


+0.488 ± 0.073 


+0.845 + 0.115 


0.144 


0.28 


FUV - NUV 


-1.312 + 0.097 


-1.835 ± 0.192 


0.237 


0.61 


-1.580 + 0.142 


-1.112 + 0.224 


0.279 


0.46 


u* - g 


+0.202 ± 0.024 


±0.373 ± 0.047 


0.058 


0.43 


+0.288 ± 0.031 


+0.307 ± 0.049 


0.062 


0.36 


9 - i 


±0.160 ±0.012 


±0.249 ± 0.023 


0.029 


0.62 


+0.225 ± 0.016 


+0.225 + 0.025 


0.031 


0.58 



its range. Using T 2 to incorporate the errors, we account for over 
60 per cent of the variance in FUV — i, FUV — NUV and g — i, 
but less than half of the variance in u* — g and NUV — i. 

The second row of panels in Figure [2] shows the equivalent 
correlations with metallicity. These are in the same sense as the age 
trends: colours which redden with increasing age also redden with 
increasing metallicity. The coefficients are significant at the ~5er 
level for all colours except FUV — i (2.6a). Again, measurement 



errors contribute substantially to the scatter, since errors on Mg/H 
individually are large. Based on T 2 , the metallicity trends and mea- 
surement errors account for about 40 per cent of the variance in 
most colours, though only 20 per cent in FUV — i. 

Above, we have modelled the colours as a function either of 
age or of metallicity, but these two parameters are not truly in- 
dependent. In particular, at fixed velocity dispersion or luminos- 
ity, older galaxies tend to be more metal poor (e.g. Trager et al. 
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2000; Smith et al. 2009). This anti-correlation appears to have 
an intrinsic component, but there is also a contribution from er- 
ror covariance, due to the 'tilt' of model grids relative to the 
line-strength indices. Regardless of its source, the mutual corre- 
lation of age and metallicity makes it difficult to determine, from 
single-parameter fits, which parameters are responsible for driv- 
ing the observed colours. Moreover, population synthesis mod- 
els show that optical colours, at least, are sensitive to both age 
and metallicity in similar degree. Hence to disentangle the effects 
of age and metallicity, we fit two-parameter models of the form 
p = ao + a t log(£ssp/10 Gyr) + az [Mg/H], to identify an im- 
proved predictor of each colour that combines the age and metal- 
licity effects to minimize the scatter. 

In the lower panels of Figure [2] we express this combination 
as log(£ssp/10Gyr) + S[Mg/H], where the ratio S = a z /a t 
differs for each colour. For the optical colours u* — g and g — i, 
we recover S = 1.6 — 1.8, in agreement with S = 1.6 expected 
from the 'three-halves rule' of Worthey (1994). Note, however, that 
Worthey's results were derived from the colours predicted by pop- 
ulation synthesis models with known age and metallicity. Here, we 
have instead derived this ratio on the basis of data for colour, age 
and metallicitfl For the UV colours, we find that NUV — i is 
more sensitive to metallicity than to age, compared to the optical 
colours (S 1 ~ 3.1), while FUV — i is somewhat more age sensitive 
(S ~ 1.0). All of the colours studied show significant dependence 
on both age and metallicity in the simultaneous fits. In Table [4] we 
compare the fit results obtained if Mg/H is replaced with Fe/H as 
the metallicity indicator. For most colours, there is little change to 
the fit results, although the scatter is generally slightly increased 
when using Fe/H. For NUV — i, however, we recover a slightly re- 
duced scatter (marginally significant), and a steeper dependence on 
age, when using Fe/H. This is likely due to strong line blanketing 
in the NUV region, which is thought to be dominated by iron lines 
(Peterson, Dorman & Rood 2001). 

Figure [3] presents these trends in an alternative manner, as 
maps of average galaxy colour in the age-metallicity plane. The 
slope of the contours in the maps reflects the relative sensitivity of 
each colour to age and metallicity. The map for NUV — i appears 
to show distorted contours indicating a more complex behaviour 
for this colour. We argue in Section l4~3l that this is due to 'contam- 
ination' of the NUV by the old hot stars responsible for the UV 
upturn. 

In most cases, the simultaneous fit to age and metallicity ac- 
counts for substantially more of the total variance in colour than 
either parameter individually, and the residual scatter is compara- 
ble to that around the colour-cr relations. The effect of measurement 
errors is smaller than for the age or metallicity taken separately, be- 
cause the age errors are anti-correlated with the metallicity errors, 
and so tend to compensate in the predicted colours. 



3.4 Systematics 

In this section we consider possible systematic error sources, and 
test whether they can make significant contributions to the residual 
scatter around the fits. 



3 The fact that our SSP parameters were derived via comparison to synthe- 
sis models does not make this a circular argument, since only spectroscopic 
indices were used to fit the ages, with no information about the broadband 
colours. 



3.4.1 Calibration uncertainties 

The aperture magnitude errors adopted for the analysis in Sec- 
tion [33] do not include contributions from systematic calibration 
uncertainties, because most galaxies in the sample are drawn from 
only a small number of observations (GALEX tiles, SDSS drift-scan 
'runs', MegaCam pointings). Hence calibration errors are likely 
to affect large sections of the dataset coherently, rather than on a 
galaxy-by-galaxy basis. In this section, we assess the evidence for 
constant shifts between individual 'observations', and their impact 
on the colour scatter, by allowing additional offset terms in the fits. 

For the g — i colour, where SDSS is the only data source, 
~90per cent of the galaxies were observed in only two drift-scan 
runs (numbers 5115 and 5087). Including a term representing off- 
sets between runs, we find that colours measured from run 5087 
are on average bluer by 0.009 ± 0.005 mag than those measured 
from run 5 1 15, at similar age and Mg/H. This is consistent with the 
claimed calibration errors of 0.01 mag in g and i (Padmanabhan et 
al. 2008). Although this term is marginally significant, its inclu- 
sion does not reduce the residual scatter around the g — i fit, which 
remains 0.029 mag. A similar test allowing offsets at the level of 
SDSS camera columns within each run does not reveal any further 
significant systematic effects at above the ~0.01 mag level. 

To test for calibration shifts in the UV versus optical colours, 
we fit a similar offset term for each GALEX tile. Here ~80per 
cent of the galaxies are drawn from two tiles (GI5_025001_COMA 
and GI2_046001_COMA3). The offset terms are not significant: 
GI2.046001_COMA3 yields colours on average bluer by 0.04 ± 
0.06 mag in FUV - i and redder by 0.01 ± 0.03 mag in NUV - i, 
for galaxies of similar age and metallicity. Again, these limits are 
consistent with the nominal GALEX calibration uncertainties of 
0.05 mag (FUV) and 0.03 mag (NUV) quoted by Morrissey et al. 
(2007). Allowing offset terms does not reduce the residual scatter. 

Applying a similar approach to the u* — g colours reveals that 
significant systematic errors are present in the MegaCam photom- 
etry. Allowing offset terms associated with each of the nine Mega- 
Cam pointings, we find that the pointing covering the south-west 
of Coma (contributing 12 per cent of the sample galaxies) yields 
colours 0.09 ±0.01 bluer in u* — g than the central pointing (57 per 
cent of the sample), for galaxies of similar age and metallicity. That 
this is due to calibration errors in MegaCam, rather than SDSS, is 
confirmed by null results obtained fitting SDSS run offsets (as for 
g — i), either separately, or in combination with MegaCam pointing 
offsets. The possibility that the offset in the south-west MegaCam 
pointing is 'astrophysical' rather than systematic can be rejected 
given that no corresponding shifts are seen for south-west region 
galaxies in residuals for the other colours (g — i, NUV — i, etc). 
Allowing for MegaCam offsets reduces the residual scatter in u* —g 
from 0.058 mag to 0.049 mag and the explained fraction of variance 
increases to R 2 = 0.63 (from 0.44). 



3.4.2 K-corrections 

We did not apply fc-corrections to the measured colours, due to the 
uncertainty in the SEDs, especially in the UV. Instead, for each 
colour we have tested for contributions due to band-shifting, by 
explicitly including an extra model term proportional to redshift. 

A significant effect (~5(J level) is recovered only for the u* — 
g colour. The redshift dependence of u* — g is not unexpected, 
since the red cut-off of the it* filter is close to the 4000 A break 
shifted to the mean recession velocity of Coma. Allowing for this 
dependence increases the explained fraction of variance in u* — g 
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by a further 6 per cent. Including both the redshift term and the 
MegaCam systematics discussed above, the explained fraction of 
variance reaches R 2 = 0.69, even without allowing for random 
errors. 

While marginal (2-3a) fc-correction terms are also detected in 
the residuals for g — i and FUV — NUV , they contribute negligibly 
to the variance (1-3 per cent) in these colours. 

3.4.3 Aperture mismatch 

Our fits compare colours derived in lOarcsec diameter photomet- 
ric apertures against age and metallicity estimates derived from 
spectroscopy in much smaller apertures (3 arcsec diameter in most 
cases). An obvious concern is whether mismatch between these 
apertures contributes to the scatter. This is a particular concern in 
the UV, which would be sensitive to any ongoing star-formation 
occurring beyond the spectrograph fibre (e.g. Salim & Rich 2010). 
As for photometric calibration shifts and fc-corrections, we can con- 
strain these effects by extending our fit with additional terms relat- 
ing to aperture mismatch. 

Since the GALEX resolution does not allow use of smaller 
apertures, we use the SDSS photometry to define a colour shift 

~ i.9 ~ *)io" — (j9~ *)s"' O n avera g e > A g i is negative, indicat- 
ing that the sample galaxies are typically bluer in g — i within the 
photometric aperture than within the spectroscopic aperture. This 
is the sense expected from metallicity gradients in purely passive 
galaxies. 

We constrain the effect of aperture mismatch in our fits by in- 
cluding an additional term proportional to A g i. For FUV — i, the 
recovered coefficient is not statistically significant, at —3.1 ± 2.1. 
The negative coefficient indicates that galaxies which are bluer in 
the photometric aperture than in the fibre in g — i are globally redder 
in FUV — i than average for galaxies of similar age and metallicity. 
This is opposite to the sense expected if the FUV — i colour scatter 
were dominated by star formation beyond the spectroscopic fibre. 
Including the A g i term in fit has negligible impact on the residual 
scatter which remains 0.29 mag. For NUV — i, the recovered coef- 
ficient is positive but not significant (+1.4 ± 1.1), and its inclusion 
does not reduce the scatter around the fit, which remains 0.15 mag. 
For the optical colours, there is a small positive correlation with 
A g i, with no associated reduction in scatter, similar to the case for 
NUV — i. Repeating these tests using colour shift in u* — g (i.e. 
A U g instead of A g i) as an indicator of aperture mismatch yields no 
further evidence for star formation beyond the fibre aperture, nor 
any non-negligible reduction in scatter around the fits. 

We conclude that unidentified star formation beyond the spec- 
trograph fibre contributes negligibly to the observed colour scatter 
at fixed age and metallicity. 

3.4.4 Internal extinction 

Although we cannot exclude the possibility that variations in inter- 
nal extinction impose some additional scatter in the relationships, it 
is clear that this does not dominate the residuals. For example, the 
residual scatter of 0. 15 mag in NUV — i, if attributed wholly to ex- 
tinction, would imply a minimum rms of 0.038 mag in g — i, which 
is far in excess of the 0.029 mag observed residual scattefl. More- 
over, the residuals in FUV — i are anti-correlated with those in 

4 This calculation assumes the Cardelli, Clayton & Mathis (1989) Ry = 
3.1 extinction law. 



g — i: galaxies that are redder than average (for their age and metal- 
licity) in the optical are unexpectedly blue in FUV — i. Such a rela- 
tionship cannot be generated only by internal extinction variations. 
Finally, the distribution of the UV versus optical colour residuals 
is skewed, with a tail of galaxies extending towards bluer colours. 
This could be caused only by these galaxies having unusually low 
dust content relative to the average, while the dispersion of internal 
extinction about its average value remains small. We consider this 
unlikely. 



3.4.5 Summary 

In summary, having reviewed the likely sources of systematic er- 
rors, we find evidence for substantial effects only in the it* — g 
colours, due to photometric offsets and band-shifting with redshift. 
Including these systematic errors, the correlations with spectro- 
scopic age and metallicity account for 70-80 per cent of the vari- 
ance in the optical colours and in FUV — NUV, but still only 
40-60 per cent of the variance in FUV — i and NUV — i. 



3.5 Residual correlations with abundance ratios 

We have shown that, using the measured ages and metallicities, 
we are able to predict optical/UV colours of passive galaxies with 
the same level of precision as using the velocity dispersion. This 
suggests that the trends of mean colour with a are adequately re- 
produced by the trends of mean spectroscopic age and metallicity 
with a. However, for FUV — i and NUV — i, where much of the 
scatter in colour was unexplained by correlation with velocity dis- 
persion, this scatter was not reduced by using models based on the 
spectroscopic constraints. 

Since element abundances in red-sequence galaxies cannot be 
adequately described by a single overall metallicity (e.g. Worthey 
1998), the abundance ratios of individual elements are obvious can- 
didates for additional parameters which might contribute to driving 
the colours. In the case of the UV upturn, in particular, the correla- 
tion of F£/V-vs-optical colours with magnesium line strength in- 
dices, and relative weakness of the correlation with iron-dominated 
indices, has led many authors to speculate that the UV upturn de- 
pends primarily on a-element abundance (O'Connell et al. 1999; 
Carter et al. 2011; Bll). 

We have tested for correlations with the element abundances 
via linear fits of the colour residuals (from the age-metallicity mod- 
els) against the ratios Fe/Mg, Ca/Mg, C/Mg and N/Mg. (We express 
all the ratios relative to Mg, rather than the more usual Fe, since 
we have adopted Mg/H as our overall metallicity indicator.) We fit 
simultaneously to all the ratios, to disentangle more cleanly the ef- 
fects of individual elements. The results are summarized in Table[5] 

We find only marginal correlations of colours with abun- 
dance ratios in most cases. For the Fe/Mg ratio, we find that 
FUV - i colour reddens by 0.07 mag, and FUV - NUV red- 
dens by 0.05 mag, for a 0. 1 dex increase in Fe/Mg at fixed age and 
Mg/H (i.e. bluer UV-vs-optical colours for Mg-enhanced popula- 
tions). However, these correlations are only ~2.5<r effects. 

Among the other residual trends with abundance ratios, the 
only significant correlations are for C/Mg: a 0. 1 dex enhancement 
in this ratio (at fixed age and Mg/H) is associated with a 0.01 mag 
reddening of g — i, a 0.03 mag redenning in NUV — i and a 
0.06 mag shift bluewards in FUV - NUV (all at 3-4tr). These 
residual trends are shown in the first two panels of Figure [4] 

To test the robustness of these results, we have also used an 
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Table 5. Residual correlations with abundance ratios. Note that because these are fits to the residuals from the age-metallicity fits, the coefficients describe 
the effects of the given abundance ratio at fixed Mg/H. Coefficients that are significant at the > 3cr level are highlighted with bold type. The final columns 
show the small reduction in total scatter when the abundance ratio trends are included in the fit. 



colour residual 


[Fe/Mg] 


[C/Mg] 


[N/Mg] 


[Ca/Mg] 


original rms 


new rms 


A(FUV - i) 


+0.67 ±0.26 


-0.33 ± 0.20 


-0.12 + 0.18 


-0.15 + 0.30 


0.291 


0.284 


A(NUV - i) 


+0.22 ±0.13 


+ 0.30 + 0.10 


+0.23 + 0.09 


+0.00 + 0.15 


0.153 


0.141 


A(FUV - NUV) 


+0.46 + 0.21 


0.63 + 0.16 


-0.36 + 0.14 


-0.15 + 0.23 


0.237 


0.222 


A(u* - g) 


+0.02 + 0.05 


+0.10 + 0.04 


+0.07 + 0.04 


-0.02 + 0.06 


0.058 


0.056 


A(s - t) 


+0.00 + 0.02 


+ 0.10 + 0.02 


+0.04 + 0.02 


+0.01 + 0.03 


0.029 


0.025 



alternative approach in which individual abundance ratio terms are 
incorporated into the fit for each colour as extra parameters, in ad- 
dition to age and metallicity. Including only Fe/Mg as a third pa- 
rameter, we find significant trends for FUV — i (+0.80 ± 0.30), 
NUV - i (+0.69 ± 0. 16) and g - i (+0.11 + 0.03). Repeating this 
test with C/Mg as a third parameter instead of Fe/Mg, significant 
trends are found for NUV - i (+0.33 ± 0.09), FUV - NUV 
(-0.45 ± 0.15) and g - i (+0.10 ± 0.02). For u - g, a significant 
trend with C/Mg (+0.16 ± 0.03) is recovered after accounting for 
the systematic effects noted in Section [3~4l These results are qual- 
itatively similar to those obtained from the simultaneous fit to the 
residuals, although individual coefficients are sensitive to the fitting 
treatment adopted. 

To summarize, decoupling the effects of multiple abundance 
ratio parameters, in addition to the age and metallicity trends, re- 
mains difficult. The strongest abundance ratio effects in the UV 
seem to be related to Fe/Mg and C/Mg, with little dependence on 
Ca/Mg or N/Mg. The abundance ratio effects make modest contri- 
butions to the total observed scatter, increasing R 2 by 8 per cent for 
NU V - i and g - i, and by 4-6 per cent in FU V - i, FU V - NU V 
and u* — g. 



4 DISCUSSION 

4.1 The index-vs-colour relations reconsidered 

As reviewed in the Introduction, early work revealed a strong cor- 
relation between the FUV-vs-optical colour and the Mg2 line index 
(Burstein et al. 1988), from which it was inferred that the FUV 
output of passive galaxies depends mainly on metallicity. Equiva- 
lent conclusions were drawn by Bll, using much improved data, 
from the strong correlation of the Mgb index against FUV — V 



colour. In this section we address the apparent disagreement be- 
tween this view and the results of our analysis which instead favour 
a dominant age dependence of the FUV-vs-optical colour, and only 
a weaker residual trend with metallicity. We base our argument on 
comparison with Bll, but qualitatively the same case applies to 
Burstein et al. (1988) and other studies following a similar analysis 
approach. 

In Figure[5] we show the correlations of FUV — g (to match 
Bll's FUV — V) with the Mg6 and H/3 line strength indices. As 
before, we treat colour as the dependent variable, and fit the rela- 
tions by minimizing scatter in the colour direction. For comparison 
we show also the fits quoted by Bll. Their approach minimizes 
scatter in the line-strength direction and their fits are made only 
to galaxies with H/3 < 1.8 A. To account for the different colour 
definitions used, the Bll fits have been shifted to pass through the 
average position of our H/3 < 1.8 A galaxies, but the slopes are 
as quoted in their paper. As always for correlations with substan- 
tial scatter, the fits differ significantly according to the direction of 
minimization. Fitting our data in the same way as Bll fit theirs, 
we would recover similar slopes. Hence at the level of the data, our 
index-vs-colour correlations are consistent with earlier work. 

The origin of the different conclusions we reach, with respect 
to Bl 1, lies in three differences in our analysis: the choice of min- 
imization direction, the treatment of the high H/3 data-points and 
the interpretation of the recovered trends. 

Minimizing in the colour direction implies asking which pa- 
rameters are the best predictors of colour. When considering age 
and metallicity as predictors (as in Section [331 there are strong 
grounds for assuming a causal relationship, since these parameters 
at least in part determine the stellar content, and hence the broad- 
band colours. However, even individual line-strength indices are 
related to the properties of the optically-dominant stellar popula- 
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tion in a fairly direct way. The UV colours by contrast have long 
appeared to exhibit large scatter with respect to all of the optical 
properties. Hence colour is the quantity with 'unknown' behaviour, 
and consequently we prefer to treat it as the dependent vari- 
able. 

The effect of excluding strong-H/3 (i.e. young and/or metal- 
poor) galaxies itself depends on the minimization direction. Mini- 
mizing residuals in colour, as in our approach, the derived trend is 
unchanged by such a cut, to first order, even when fitting H/3 itself. 
By contrast, if we minimize residuals in index strength (as done by 
Bll), the resulting fit is highly biased for H/3, since all H/3 val- 
ues are included at blue colours, but only unrepresentatively low 
H/3 values contribute at redder colours. The bias is less strong for 
Mg6, since the selection is not in the fitted quantity, but still present 
because Mgb tends to be anti-correlated with H/3. (Note that our 
fits may be slightly affected by a bias in the opposite sense: since 
our sample is limited by magnitude error in the FUV, the reddest 
galaxies may be preferentially excluded at the faint end.) 

Thus we can understand why the different choice of fit di- 
rection and exclusion of strong-H/3 galaxies led Bll to a steeper 
slope for the colour-vs-Mg6 relation (inverting the values in their 
table2gives — 2.51q "7 mag A" 1 compared to our slope of —0.63 ± 
0.05 mag A -1 ) and a dramatically steeper slope for the colour- 
vs-H/3 relation (+8.3^2 5 mag A -1 compared to our +0.97 ± 
O.lOmagA -1 ). It might be expected, then, that Bll would infer 
a much steeper age dependence for the FUV-vs-optical colour than 
we do. Instead, they imply that the trend results from the known 
dependence on H/3 on metallicity, and hence that a correlation be- 
tween FUV and metallicity drives both the colour-Mgb and colour- 
H/3 relationships. Bll do acknowledge that the high-H/3 outliers 
from their fits are driven by younger ages. Our interpretation of 
their results is that these 'outliers' in fact fall along the more mean- 
ingful correlation obtained by fitting the full dataset. 

In summary, the conclusions reached in our work disagree 
with those of B 1 1 mainly through our preferred approaches for fit- 
ting and interpreting the data. 



4.2 The origin of the FUV upturn 

One motivation for this work was to explore the age and metallicity 
dependence of the FUV upturn, as a method to distinguish between 
the 'metal-rich single-star hypothesis' (e.g. Yi et al. 1997b) and the 
'binary-star hypothesis' (Han et al. 2007) for the origin of the FUV 
flux in passive galaxies. Although only ~13 per cent of galaxies in 
our sample have classic FUV upturns, in the sense of having rising 
/a-vs-A below 2000 A (FUV - NUV < 0.9), we prefer a 'con- 
tinuous' treatment of the colours, rather than imposing a discrete 
classification as done for instance by Yi et al. (201 1). 

In both the Han et al. and Yi et al. models, the FUV up- 
turn is caused by helium-burning stars with very thin hydrogen 
envelopes, having effective temperatures ~25 000 K. In addition 
to this extreme horizontal branch (EHE0) phase itself, the subse- 
quent post-HB evolution of such stars can include long-lived UV- 
luminous phases (e.g. the failed asymptotic giant branch, or AGB- 
Manque). These stars were identified as probable sources of the 
UV upturn through their good match to the spectra (both contin- 
uum and lines) of strong-upturn galaxies like NGC 1399 (e.g. Fer- 
guson et al. 1991). Hubble Space Telescope imaging of M32 by 



5 Han et al. (2007) refer to these stars more generally as hot sub-dwarfs, 
since they consider a wide range in formation mechanisms. 



I 

> 
z> 




— B1 1 fit (minimizes horizontal scatter) 
Our fit (minimizes vertical scatter) 



2.5 3.0 3.5 4.0 4.5 
Mgb5177 



5.0 



5.5 




Figure 5. The linestrength-vs-colour relations, for comparison to previous 
work, in particular Bureau et al. (2011, Bll). The green line shows our fit 
to our data, minimizing in the colour direction, as throughout our work. 
The red line shows the Bll fits (to their data), which minimize in the line- 
strength direction and exclude galaxies with H/3> 1.8 A. 



Brown et al. (2000, 2008) has confirmed the presence of EHB and 
post-EHB stars, and their dominance of the FUV emission in this 
(weak-uptum) elliptical. 

The Han et al. and Yi et al. models differ in the channels by 
which stars reach the EHB, and hence on their dependence on age 
and metallicity. In the single-star model, Yi et al. (following Greg- 
gio & Renzini 1990) postulate EHB formation through an enhanced 
efficiency of mass loss on the red giant branch at high metallic- 
ity, leading to low envelope mass and high temperature on the HB. 
The evolution of such stars after they leave the HB is also sensitive 
to envelope mass, and hence to metallicity, with a larger fraction 
of high-metallicity stars following the UV-luminous AGB-Manque 
pathway instead of evolving onto the asymptotic giant branch. Fi- 
nally, there is a strong age dependence because stars of smaller 
initial mass need to lose less mass to reach the EHB and the UV- 
bright post-HB phases. In the models favoured by Yi et al. (1997b, 
their Figure 16), the envelope-mass threshold governing the post- 
HB evolution causes an extremely rapid onset of the UV upturn, 
e.g. a 2^1 mag boost in FUV between ages of lOGyr and 15 Gyr. 
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For younger ages the evolution is much slower, perhaps equivalent 
to a slope of 0.3 — 0.6 mag per decade in age. The effect of metallic- 
ity is primarily to affect the age at which the rapid UV boost occurs, 
such that metallicity trends at constant age may be very strong or 
very weak, depending on the age range considered. 

In the Han et al. scenario, EHB stars form through a vari- 
ety of binary interactions, including stable Roche-lobe overflow, 
common-envelope evolution and mergers of binary white dwarfs. 
In these models, the FUV-vs-optical colours are only weakly de- 
pendent on age for an SSP. For example, from Table 2 of Han et 
al. (2007), fitting the FUV — r colour over a range of 3 Gyr to 
12Gyr, we obtain a slope of —0.3 mag per decade in age. The in- 
fluence of metallicity is not fully addressed: although channels for 
formation of EHB stars are not thought to depend on metallicity, 
their subsequent evolution may do so, and this is not followed in 
the model. 

In testing the observed age and metallicity dependence of 
the upturn, we note that the proposed correlations are 'indirect': 
the hot stars responsible for the FUV excess are too faint at opti- 
cal wavelengths to contribute substantially to the spectroscopic in- 
dices. Rather, we are testing the hypothesis that (a) the mechanisms 
for producing these stars are affected by the age and metallicity of 
the stellar population from which they evolved, and (b) that this 
population was either the same as, or at least correlated with, the 
stellar population which today dominates the optical spectrum. 

We adopt the observational results for FUV — i as the most 
direct tracer of the upturn population (see Section |4~4l for comments 
on FUV — NUV). The fits presented in Section [331 suggest that 
FUV — i is primarily sensitive to age, becoming bluer for older 
galaxies with a slope of 1.1 ±0.1 mag per decade in age. At fixed 
age, there is also a dependence on metallicity, such that FUV — i 
becomes bluer by 1.1 ± 0.2 mag per decade increase of Mg/H. The 
sense of the observed age trend is as predicted by both the Han 
et al. and Yi et al. models, but its slope is much steeper than pre- 
dicted by Han et al. As noted above, the age dependence in Yi et 
al. is non-linear, becoming extremely steep at ages above 10 Gyr. 
Some steepening is evident in the upper left panel of Figure[2]as an 
increased scatter on the blue side of the distribution at ages above 
~ 10 Gyr. Quantifying this by splitting the sample at the median age 
of 9 Gyr, we find that the older galaxies show a slope of 1.5 ± 0.3 
per decade in age, while the younger galaxies have only a marginal 
dependence 0.5 ± 0.2. This behaviour is difficult to reconcile with 
the Han et al. model, but is a generic prediction of the Yi et al. sce- 
nario. We note that a strong dependence at the oldest ages is also 
suggested by the apparently rapid redshift evolution of the UV up- 
turn (e.g. Brown et al. 2003; Ree et al. 2007, but see also Atlee, 
Assef & Kochanek 2009). The Yi et al. model also predicts that 
the FUV-to-optical colour becomes bluer with increasing metal- 
licity, but the strength of the metallicity effect itself depends on 
age. Based on figure 16 of Yi et al. (1997a), we would expect a 
trend of order 10 mag per decade in Mg/H for galaxies with ages of 
15 Gyr, but only ~0.8 mag per decade at 5-10 Gyr ages. The latter 
is broadly consistent with the trend derived from our fits. 

Despite recovering a strong age dependence in the FUV — i 
colour, we continue to find a large scatter around the average trend, 
with variations in the optically-dominant stellar populations ac- 
counting for only -~39 per cent of the total variance in colour, and a 
further 19 per cent is contributed by measurement error. There are 
only marginal correlations of FUV — i with element abundance ra- 
tios, which do not substantially reduce the scatter. Hence a simple, 
deterministic, explanation of the FUV colours appears inconsistent 
with SSP models, at least to the extent that they can be constrained 



from the optical spectra. Considering more realistic extended star- 
formation histories (SFH), in combination with the threshold age 
effect in the Yi et al. models, perhaps offers a more promising ex- 
planation for the scatter in FUV — i. The SSP-equivalent age is 
essentially a weighted average over the true distribution of stel- 
lar ages. Because in the Yi et al. scenario the oldest populations 
produce dramatically more FUV flux, the FUV — i colours for a 
galaxy with extended SFH will be sensitive to the fraction of stars 
in the old tail of the age distribution. As a simplified example, con- 
sider a case where populations older than 12 Gyr receive an instan- 
taneous 2 mag boost to their FUV flux. A galaxy with gaussian SFH 
with mean age to=lO Gyr and dispersion ot=\ Gyr has 2.3 per cent 
of its stellar content FUV-boosted, resulting in a total flux 0.12 mag 
brighter in FUV than an SSP of the same age. Keeping the same 
mean age, but increasing the dispersion to <7 t =1.5Gyr results in 
boosting ~9.1 per cent of the stellar content, increasing the total 
flux by a further 0.31 mafl Thus if the FUV luminosity increase 
for the oldest galaxies is sufficiently strong and rapid, then the ob- 
served excess scatter of ~0.3 in FUV — i could be attributable 
to a modest variation in star-formation time-scale, at a given mean 
ag4H Further investigation is required to test whether this explana- 
tion is viable given more realistic and detailed models of the flux 
evolution. 

In conclusion, the observed relationship between FUV — i 
colour and spectroscopic age and metallicity appears to favour the 
metal-rich single-star EHB hypothesis for the origin of the UV up- 
turn. This scenario generically predicts a dependence on age in par- 
ticular among the oldest galaxies, and on metallicity. The large scat- 
ter at given SSP-equivalent age and metallicity can possibly be ex- 
plained through galaxy-to-galaxy variations in early star-formation 
history. The binary hypothesis cannot account naturally for these 
results. 



4.3 Excess scatter at NUV: residual star-formation, or 
UV-upturn leakage? 

A second motivation for our study was to investigate the cause of 
the large scatter seen in NUV-vs-optical colours. 

The NUV-to-optical scatter has been interpreted as evidence 
for recent or widespread ongoing star formation in optically-red 
galaxies (e.g. Kaviraj et al. 2007). Sample definition is critically 
important to this result, since objects with ongoing star formation 
will evidently have blue UV colours. For our analysis, we excluded 
galaxies with central Ha emission from the spectroscopy. The re- 
maining objects have only ~0. 16 mag scatter around the NUV — i 
colour-magnitude (or colour— a) relation, compared to values in 
the range 0.3-0.6 mag reported in the literature (Yi et al. 2005; 
Boselli et al. 2005; Haines, Gargiulo & Merluzzi 2008; Rawle et 
al. 2008). Star formation occurring further from the nucleus can- 
not be excluded, but does not appear to dominate the colour scatter 
(Section |3A3l 

We now consider the correlations of NUV — i with spectro- 
scopic age and metallicity. Unlike in the FUV case we are here test- 

6 In general, for a gaussian SFH with mean age to, age dispersion at, and 
a flux-boost factor b imposed on ages above T, the total flux is increased by 
a factor /= | (1 + 6) + |(1 - 6) erf (^M). 

7 However, if the a element abundance ratio is interpreted as a star- 
formation time-scale indicator (Thomas et al. 2005), then this explanation 
also predicts that bluer FUV — i (at given age) would be accompanied by 
higher Fe/Mg ratios, contrary to the observed trend. 
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ing 'direct' correlations, in the sense that the same stars (in partic- 
ular warm main-sequence turn-off stars) are expected to contribute 
to the NUV colour variations and to the spectroscopic ages. We find 
significant trends for bluer colours at lower metallicity (the domi- 
nant effect, as found also by Rawle et al. 2008) and at younger 
age. While the age dependence is qualitatively as expected from 
evolutionary synthesis models, the slopes that we obtain from the 
fits (e.g. ~0.2mag per decade in age) are much shallower than 
the predictions derived by Dorman et al. (2003), who report a 2- 
3 mag change in NUV — V per decade in age. The evolution of 
NUV — i in the Bruzual & Chariot (2003) models is more grad- 
ual, with ~0.65 mag change per decade in age, though still steeper 
than our measured trends. (The Bruzual & Chariot models do not 
include blue horizontal branch stars; their presence would further 
dilute the predicted age trend.) 

Age and metallicity variations account for only 20 per cent 
of the total variation in NUV — i, leaving a residual scatter of 
0.15 mag around the fit. Splitting the sample by spectroscopic age, 
we find that the fit is poorer (i.e. larger scatter and smaller R 2 ) 
for the oldest galaxies in the sample. For galaxies older than the 
median age of 9Gyr, the age and metallicity effects account for 
only 12 per cent of the variance. For galaxies younger than 9Gyr, 
correlations with age and Mg/H jointly account for 44 per cent of 
the total variance. 

As noted elsewhere, the SSP-equivalent is only a weighted 
average over the unknown true SFH of the galaxy. For example, 
an SSP age of ~5 Gyr may represent a population formed mostly 
at early epochs, but with a small mass-fraction added in the past 
~Gyr. We have tested whether the remaining scatter in NUV — i 
is due to young sub-populations, using the HSf index. Discrepan- 
cies in HSf compared to the (H/3-based) SSP-equivalent age es- 
timates can be interpreted as evidence for star-formation within 
the past ~2.5 Gyr (Serra & Trager 2007), or alternatively as ev- 
idence for contributions from other warm stars, e.g. on the blue 
HB (Schiavon et al. 2004). We consider the residuals from the fit 
of NUV — i versus age and metallicity, and compare them to the 
residuals from a fit of HSf against age and metallicity (third panel 
of Figure [4] The residuals from the two fits are anti-correlated, 
as expected: galaxies with stronger HSf than expected from their 
SSP-equivalent ages and metallicity are also bluer than the average 
galaxy at that age and metallicity. The correlation is significant at 
the ~4<r level, and this dependence accounts for ~8per cent of the 
total variance in NUV — i. Thus, variations in true star-formation 
history (or blue HB content), as traced by HSf make a measurable, 
but minor, contribution to the remaining NUV — i scatter at a given 
SSP-equivalent age and metallicity. (Equivalent tests for the other 
colours reveal 2-3<r trends for u* — g and g — i, which < 3 per cent 
to the variance in these colours.) 

The above results suggest that the excess scatter at NUV is 
related to old populations, rather than to young stars. We hence ex- 
plore here whether NUV — i is significantly affected by the old 
helium-burning stars, i.e. the UV-upturn sources. The possibility 
that such sources make a significant contribution to the flux at NUV 
has previously been discussed in previous work (e.g. Dorman et al. 
2003; Rawle et al. 2008). If we adopt FUV — i as an indicator for 
the strength of the upturn population, we find that A(NUV — i), 
defined as the residual from the age-metallicity fit for this colour, is 
correlated with FUV — i at the 6<r level (fourth panel of Figure[4]l. 
(By contrast the NUV — i colour itself shows only a ~2a correla- 
tion with FUV — i, because other factors are important in driving 
the NUV flux.) This result suggests that some of the excess scatter 
in the fits to NUV — i results from 'leakage' of the FUV upturn. 



Incorporating FUV — i as an additional term in a simultaneous 
fit for NUV — i reduces the scatter from 0.15 mag to 0.13 mag. 
In the resulting fit, the coefficient of FUV - i is 0.31 ± 0.04; the 
slope with log Tsgp is increased to ~0.55 (in better agreement with 
the Bruzual & Chariot models), and that of [Mg/H] is increased to 
~1.0. The resulting coefficient of determination is increased from 
R 2 = 0.20 to R 2 = 0.48. The effect of 'correcting' NUV - i 
for the FUV leakage can be seen in the upper-right panel of Fig- 
ure^ Comparison to the equivalent panel for uncorrected NUV — i 
shows that the correction makes the oldest galaxies redder (since 
these are the galaxies with strongest UV upturn) and hence yields 
a map similar in form to the optical colours. It has been suggested 
that the FUV-upturn sources may even affect the 'optical' U band 
(e.g. Yi et al. 1997b). However, we find no correlation of the u* —g 
(or u* — i) residuals with FUV — i, equivalent to the trends seen 
at NUV. We infer that "leakage" of the upturn does not extend sig- 
nificantly to the fairly red u* bandpass used in our MegaCam ob- 
servations. 

Thus, whatever the sources of the FUV variation might be, 
they are likely to cause a significant fraction of the scatter observed 
in the NUV as well. Of course, this result does not, on its own, 
prove that the correlated excess variation in FUV and NUV are both 
due to old hot stars (i.e. the classic UV upturn population), rather 
than both being due to residual ongoing star formation. We consider 
the latter possibility unlikely though, because on other grounds it 
is clear that the UV upturn itself is in general not due to young 
stars (O'Connell 1999, and references therein). Moreover, recall 
that FUV — i is bluer for older spectroscopic age. If ongoing star 
formation were responsible for variations in FUV — i, it would 
have to occur preferentially in galaxies that have not formed stars 
at the intermediate ages (~1 Gyr) which strongly affect the SSP 
fits. 

Finally we note that Salim & Rich (2010) have presented Hub- 
ble Space Telescope FUV imaging revealing clear signatures of 
star-formation in a sample of 'quiescent early-type galaxies' that 
have no detectable Ha in SDSSjj. The Salim & Rich sample how- 
ever was constructed explicitly to select galaxies that are unusually 
blue in FUV — i, relative to their other properties, so that by defi- 
nition they are outliers from the overall population. Indeed, almost 
all of their sample galaxies are much bluer, in both FUV — i and 
NUV — i, than any of our Coma galaxies (measured relative to 
the peak of the red sequence in the respective studies). Although 
their result confirms that nuclear Ha measurements alone cannot 
exclude all star-forming galaxies, it does not imply that such galax- 
ies are common among samples selected to be representative of the 
optical red/passive sequence. 

We conclude that a majority of the observed scatter in the 
NUV-vs-optical colour can be adequately accounted for by a com- 
bination of (i) variation in age and metallicity as inferred from op- 
tical spectroscopy (~20 per cent of total variance), (ii) a signifi- 
cant contribution from the same old hot stars which dominate the 
FUV scatter (~30 per cent), (iii) a minor contributions from abun- 
dance ratio effects and variation in SFH at given SSP-equivalent 
age (~15per cent) and (iv) measurement error (~15 per cent). 
There is little need to invoke widespread 'residual star formation' 
to explain the NUV colour scatter in our sample. 



8 Although in fact half of the examples they show are evidently spirals on 
the basis of the SDSS imaging alone. 
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Table 6. Summary of the sources of variance in the colours as determined in this paper. The contributions are expressed as a fraction (in per cent) of the 
total variance in each colour. We distinguish 'astrophysical' sources of variance, i.e. correlations with intrinsic galaxy properties, from those due to systematic 
errors in the measurements. The sum of these constitute the total 'structural' component of the variance, i.e. the fraction attributable to known dependencies. 
Adding the variance due to random errors, we arrive at the total fraction of variance that has been accounted for. All components below the ~4per cent level 
are omitted for clarity. 



Contribution 



FUV ■ 



NUV - i 



FUV - NUV 



9 ~ « 



Section 



Astrophysical: 



Systematics: 



Age/metallicity 
Abundance ratios 
Young sub-populations 
UV-upturn leakage 

Photometric calib 
fc-correction 



39 

5 



20 

8 
8 

28 



62 

4 



44 



19 



63 



3.3 
3.5 
4.3 
4.3 

3.4.1 
3.4.2 



Total structural component 

Random errors 

Total explained component 



44 

19 
63 



64 

17 
81 



66 

21 
87 



75 
9 
83 



72 
9 
81 
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4.4 A comment on the tight correlations for FUV-NUV 

We have shown that the overall strength of the UV upturn, as traced 
by FUV — i, shows scatter that cannot be accounted for by vari- 
ations in the optically-dominant stellar populations. However, the 
FUV — NUV colour is surprisingly well behaved, in the sense 
that most of its variation is predictable from the velocity dispersion, 
or from age and metallicity, with a precision similar to that found 
for the optical colours. (This statement is in fact equivalent to our 
conclusion that leakage from the UV upturn contributes much of 
the excess scatter in NUV — i.) Age and metallicity variations ac- 
count for ~60 per cent of the total variance in FUV — NUV, with 
a further ~20 per cent attributable to measurement errors. 

That the FUV — NUV colour follows tighter relations than 
the FUV-vs-optical colours has been noted before (Donas et al. 
2007), and Bll conclude that FUV - NUV should be the pre- 
ferred colour to investigate the UV upturn. We disagree with this 
position: given the consensus that the upturn sources are hot old 
stars, the crucial quantity is their incidence relative to cool stars 
of similar age, which is best traced using the red-optical or near- 
infrared flux. 

Simplistically, the FUV — NUV colour should be sensitive 
instead to the temperature of the upturn sources, which would pro- 
vide further constraints on the origin of the upturn. For instance, 
Figure 14 of Yi et al. (1997b) shows that within their model, the 
post-EHB stars generate more FUV flux than the EHB itself, while 
in the NUV the contributions of EHB and post-EHB are compara- 
ble. In practice however, any attempt to extract information from 
the FUV — NUV colour alone will be confused by its sensitivity 
to at least three influences: (i) the age and metallicity of the pre-HB 
stellar population, via the main-sequence turn-off, (ii) the overall 
strength of the UV-upturn population, and (iii) the spectral slope of 
the upturn. In particular the strong observed FUV — NUV cor- 
relations likely arise from a compounding of effects (i) and (ii): 
increasing age and metallicity have the effect of suppressing the 
NUV flux from the main sequence, as well as boosting the FUV 
flux from the helium-burning populations. 

We conclude that the FUV - NUV colour will be of lim- 
ited use in deciphering the origin of the upturn, unless additional 
information from longer wavelengths is incorporated. 



5 CONCLUSIONS 

We have analysed the UV and optical colours of 150 galaxies in 
the Coma cluster, selected to lie on the (optical) red-sequence and 
to have no detectable Ho emission. Using ages and metallicities 
derived from optical spectroscopy, we have performed a purely em- 
pirical test for the origin of scatter in the colours. Our primary con- 
clusions are as follows: 

(i) All of the UV and optical colours we investigate show strong 
(> 5ct) correlations with luminosity and with velocity dispersion. 

(ii) The average trends in galaxy colours as a function of ve- 
locity dispersion can be accounted for by their correlation with the 
spectroscopically-measured SSP-equivalent ages and metallicities. 
NUV — i,u* — g and g — i become redder with increasing age or 
metallicity, while FUV — NUV and FUV — i become bluer. 

(iii) For the optical colours, and for FUV — NUV, most (70-80 
per cent) of the variance in galaxy colours can be accounted for by 
correlations with spectroscopic age and metallicity, once random 
and systematic error sources are included. 

(iv) For the FUV — i colour, which traces the UV upturn, a 
strong anti-correlation with age is observed, combined with an anti- 
correlation with metallicity. The correlations are in the sense pre- 
dicted by single-star EHB models for the origin of the FUV flux in 
passive galaxies (Yi et al. 1997b), and contrary to the weak depen- 
dences predicted by the binary hypothesis of Han et al. (2007). 

(v) However, only ^60 per cent of the variance can be ac- 
counted for, even when measurement errors are included. Moreover 
the galaxies with the oldest SSP-equivalent ages (> 10 Gyr) appar- 
ently show an excess flux in FUV which is not captured by the lin- 
ear fits. Such a steepening at high age is predicted by the single-star 
EHB models. The excess scatter at FUV could result from a vary- 
ing contribution from such ancient stars in galaxies with younger 
SSP-ages but extended star-formation histories. 

(vi) Excess scatter is also observed for the colour NUV — i. 
The residuals are highly correlated with the FUV — i colour, sug- 
gesting that the stars responsible for the UV-upturn also contribute 
significantly to the scatter in NUV. Variation in the true extended 
star-formation histories, at a given value of the SSP-equivalent age, 
makes only a small contribution to the NUV — i colour scatter. 
There is no need to invoke widespread 'residual star formation' to 
account for the observed NUV scatter in our sample. 
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(vii) At fixed age and Mg/H, colours are weakly correlated with 
the Fe/Mg and C/Mg abundance ratios, with no significant depen- 
dence on Ca/Mg or N/Mg. Abundance ratio variations, at fixed 
Mg/H, contribute only modestly to the scatter in the UV colours. 

(viii) Combining all effects addressed in this paper, we can ac- 
count for 80-90 per cent of the scatter in all of the colours studied, 
except for FUV — i. The contributions are summarised in Table|6] 

Because we have limited our analysis to empirical correla- 
tions, the above results are independent of any specific population 
synthesis models for the UV colours. The measured correlations 
and variances should provide a good comparison dataset for future 
model-based investigations. 
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